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Abstract. Polya showed in his 1921 paper that the generating function of the return probability 
for a two-dimensional random walk can be written in terms of an e|liMic integral. In this paper we 
present a similar expression for a one-dimensional quantum walk. || || 1 1 

1 Introduction 

The classical random walk (RW) is one of the most popular models to analyze many prob- 
lems in various fields. The quantum walk (QW) is expected as a counterpart of the RW in 
quantum systems. The intensive study on QWs has a short history, about ten years, how- 
ever the attention to the QW increases in the scientific community, for example, physicists, 
mathematicians, and computer scientists. Recently a book and lecture notes on the QW 
were published [H |2]. For excellent reviews, see |3l H]. There are two types of QWs, i.e., 
discrete- and continuous-time walks. The relation between them was investigated in O |6]. 
Ambainis et al. [7] intensively studied the one- dimensional (ID) discrete-time case by the 
Fourier analysis and path counting method. Recently it was shown by |H] that how the 
theory of CMV matrices gives a natural tool to study the ID discrete-time QW. Concerning 
the CMV matrix, see [HI ED], for instance. In this paper, we consider the ID discrete-time 
QW which is determined by a 2 x 2 unitary matrix. We prove that the generating function 
of the return probability for the ID QW can be written in terms of an elliptic integral by 
a path counting method. A similar expression for a 2D RW was presented by Polya in 
1921 yjj. The return probability of a 3D RW is also given by an elliptic integral and was 
evaluated by Watson in 1939 [12]. It is known that elliptic integrals and their inverses, i.e., 
elliptic functions, come from many branches of mathematics and physics. Our statement on 
the generating function and the elliptic integral may be the tip of an iceberg. Pursuing the 
study along this line would be one of the interesting future problems on the QW. 

The reminder of the paper is organized as follows. Section 2 treats the definition of the 
QW. In Sect. 3, we present our main result (Theorem 13. 2p . Section 4 is devoted to proof of 
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Proposition 13.11 Concerning the analysis on the QW, there are some methods, that is, the 
Fourier analysis, the path counting, the CMV matrix-based method. In our proof, we use a 
path counting method as in the classical 2D RW case stated in Sect. 3. Finally, we give a 
conclusion in Sect. 5. 



2 Model 

In this section we define the discrete-time QW on Z considered here, where Z is the set of 
integers. In general, the time evolution of the QW is determined by a 2 x 2 unitary matrix: 



U 



a b 
c d 



with a,b,c,d & C and C is the set of complex numbers. The QW defined by the Hadamard 
gate U = H with a = b = c = —d = 1/a/2 is often called the Hadamard walk and has been 
extensively investigated in the study of the QW: 
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In the present paper, we focus on the Hadamard walk. The discrete-time QW is a quantum 
version of the random walk with additional degree of freedom called chirality. The chirality 
takes values left and right, and it means the direction of the motion of the walker. At each 
time step, if the walker has the left chirality, it moves one step to the left, and if it has the 
right chirality, it moves one step to the right. Let define 



1^) 



1 




\R) 
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where L and R refer to the left and right chirality state, respectively. To define the dynamics 
of our model, we divide U into two matrices: 



P 



a b 




, Q 





c d 



with U = P + Q. The important point is that P (resp. Q) represents that the walker moves 
to the left (resp. right) at position x at each time step in the following. We let S„(/,m) 
denote the sum of all paths starting from the origin in the trajectory consisting of / steps 
left and m steps right at time n with I + m = n. For example, we have ^2(1, 1) = QP + PQ 
and 



^4(2, 2) = g2p2 ^ p2Q2 ^ QpQp ^ pQpQ ^ pQ2p ^ Qp2Q 



(2.1) 



In this paper, we take (p^ = "^[I/a/2, i/\/2] as the initial qubit state, where T is the transposed 
operator. Then the probability distribution of the walk starting from ip^, at the origin is 
symmetric. The probability that our quantum walker is in position x at time n starting 
from the origin with ^9^, is defined by 



P{Xn = x) = ||H„(/,m)v9*| 
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where n = I + m and x = —I + m. Then the return probabihty at the origin at time n is 
given by 

p„(0) = PiXn = 0). 

By definition, p2n+i(0) = for n > and it is enough to study only even times 2n. 

3 Result 

In this section we present our results. By definition of the Hadamard walk, we can directly 
compute 

Po(0) = l, p2(0) = ^ = 0.5, p4(0)=p6(0) = ^ = 0.125, 

g 25 
P8(0) = Pio(O) = Y^ = 0.07031 . . . , puiO) = puiO) = ^ = 0.04882 . . . , 

1225 

Pi6(0) = Pi8(0) = = 0.03738 .... (3.2) 

In general, P2n(0) can be written by Legendre polynomials, Pn{x), as follows. As for special 
functions, see [TBI. 



Proposition 3.1 po(0) = 1 and 

P2n{0) = l [{Pn-im" + {PniO)}'] {n>l). 

The proof appears in the next section. By Proposition 13. H P2n+i(0) = and 

P2n(0) = — ) , (3.3) 



n 



we get 



1 1 /2m\ ^ 

pUO) = Pm(0) = -{P2mm' = f ^ j (m > 1). (3.4) 



So (13.21) can also be obtained by (13. 4p . From now on we derive our main result (Theorem 13. 2p 
from Proposition 13.11 We begin with 



n=0 



J]p„(0)z" = 5^P2n(0)^^" = - 
n=0 n=0 

Therefore by (13. 3 p we see that 

n=0 n=0 n=0 \ / \ / 
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where K{k) is the complete elhptic integral (see [13]), i.e., 

K{k) = / , = / , < A: < 1 . 

Then the generating function of the return probability for the ID Hadamard walk can be 
expressed by K{k): 

Theorem 3.2 

°° 1 -U ^2 1 

^ — ' TT Z 

n=0 

Here we consider the generating function for the classical case. The dimensional classical 
(simple) RW whose transition probability from the origin to x is given by l/2d for any 
X = (xi, . . . , Xd) € TJ^ with Yl'k=i \^k\ = 1- Let ''^''(0) denote the return probability for the 
RW starting from the origin at time n. In the ID RW case, it is well known that 

P£^^(0) = ^(^J'), p£i(0) = (n>0). 

So we get 

CO -I 



n=0 



For the 2D RW, we find that 



p£'Ho) = {pr/^(o)}' = ^ (^;) , pa\(o) = (n > 0). 

Concerning the derivation, see for instance [14J. Therefore we have 

oo 

j]p(r'^)(o)z" = -ir(^). 

n=0 

The result was obtained by Polya in 1921 ([11], p. 160). In addition, the above expression 
is very similar to that for the ID Hadamard walk (see Theorem 13.21) . The probability that 
the 3D RW returns to its starting point, F = 1-0'^ = 0.34053 . . ., is given by K{k) in the 
following ([I5|, p. 103): 

G = 4 / k( — -I de. 



TT^ J _^ \3 — COS 6* 

We should remark that the definition of the complete elliptic integral in [15] is (2/7r) x K{k) 
in our notation. This integral was first evaluated by Watson in 1939 [T2] : 



G = 3(18 + 12^2 - lOVs - 7Vq) {K{2V3 + Vq - 2^2 - 3)}^ x {2/tt 



1.51638.... 
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Table 1: Products of P, Q, R, S. For example, PQ = bR. 





P 


Q 


R 


S 


p 


aP 


bR 


aR 


bP 


Q 


cS 


dQ 


cQ 


dS 


R 


cP 


dR 


cR 


dP 


S 


aS 


bQ 


aQ 


bS 



4 Proof of Proposition 13.1 



In this section we will prove Proposition 13. II bv a path counting method. First we introduce 
the useful matrices to compute S„(/,m): 



R 



c d 




, S 





a b 



where a = b = c = —d = l/\/2. In general, products of the matrices P, Q, R and S are 
given in Tabled] From this table and fl2.ip . we obtain 

^4(2, 2) = bcdP + abcQ + b{ad + bc)R + c{ad + bc)S. 

We should note that P,Q,R and S form an orthonormal basis of the vector space of complex 
2x2 matrices with respect to the trace inner product {A\B) = tT{A*B), where * means the 
adjoint operator. So S„(/,m) has the following form: 

m) = pn{l, m)P + qn{l, m)Q + r„(/, m)R + m)S. 

The exphcit forms of Pnil, Tn), qn{l, m),rn{l, m) and Snil, m) can be computed as follows (see 
for instance [16] flTj): 



Lemma 4.1 When I A m(:= min{Z, m}) > 1, we obtain 

n-l (;-l)Am 



qn{l,m) 



V2 

1 

7i 



7=1 

1 lA{m-l) 

E (-1 

7=1 

n— 1 lAm 



m— 7 



771—7- 



/ — 1 \ /m — 1 



7 



7-1 



^ f I — l\ f m — 1 
7-1 



rnil,m) = Snil,m) 



By using Lemma 14. 1^ we have 



1 

7! 



^(_iy«-7 
7=1 



7 

/ - 1 \ /m - 1 
7-1; V7-1 



-2n 



V2) 



2n 



(-irE(-i) 



7=1 



1 \ 



2n+l 



n — 1 
7-1 



1)^ - 1 



n n 

7 7 

+ 2 

7 



- 2 " 


1 


' 1 " 


n 
7 


7! 


i 



\ 7 V7 - 1 



n + (n — 27)^ 
— (n — 27) + m 
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Since J92n(0) = | |H2„(n, n)v9*| p, we get 



P2n(0) 



2n 



2n 



7 V7-1 



- 27) 



2 J i^'^V f n — 1 

'':^^^V7-i 

n n 

-""EE 

7=1 (5=1 

n n 



4)7+5 _ 1\ 2 _ 1 
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V7-V 



Furthermore we will rewrite P2n(0) by using the Jacobi polynomial, P^'^{x), which is orthog- 
onal on [—1, 1] with respect to (1 — xY{l + xY with z/, > —1. Then the following relation 
holds: 



''^'^^ = r(n + i)r(^+\) ^ + ^ + + + ^)/2), 



(4.5) 



where T{z) is the gamma function. Therefore we see that 



E 

7=1 



1)^-1 (n - 1 
7-1 
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2Fi(-(n-l),-(n-l);2;-l) 
2"-i2Fi(-(n-l),n + l;2;l/2) 

on— 1 

2 -^(1,0) 



n 



(4.6) 



The first equality is given by the definition of the hypergeometric series. The second equality 
comes from the following relation: 

2Fi(a, 6; c; z) = {I - z^^Fiia, c-b;c; z/{z - 1)). 

The last equality follows from (14. 5p . In a similar way, we find that 



-^(0). 



(4.7) 
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By using (14.61) and (14. 7p . we get 

P2n(0) 



2n^ X 



n 

o2(ra-l) 

f_ :p(i 



1 
2 

— 4n X 

^"-^Wf ^pa;)(o)p(»,;)(o) + |pn(o)^' 



n 



e;'(0)pr-:;'(0) + 4x^<"-') {Fi°:;'(0)}- 



= 5[{/'.;i:;'(o)-PH'(o)}V{pi°:?(o); 

From (6.4.20) in [I3]: 

(n + a + l)Pi-^'\x) -in + l)Pi:f (x) = (^!i±^±A±Ml^p(«-w./3)(3,), 
we have 

i'rHo)-pf^?(o) = prHo). 

Combining fH^ with l^i^ yields 

P2„(o) = {pr^(o)}V{pi°:°)(o)}'. 

Since Pn(0) = Pn°'°''(0), the proof of Proposition 13.11 is complete. 



(4.8) 



(4.9) 



5 Conclusion 

In the present paper, we showed that the generating function of the return probability for 
the ID Hadamard walk can be written in terms of an elliptic integral. Our expression 
corresponds to the result for the classical 2D RW given by Polya in 1921. Indeed, his 
expression is also written in terms of an elliptic integral. One of the future interesting 
problems is to find similar expressions for general ID and higher dimensional QWs. Resent 
study on the relation between discrete- and continuous-time QWs, e.g. [5l|6], would allow us 
to give a new aspect between our result and the corresponding one for the continuous-time 
case. 
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